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INTRODUCTION 


I  A  ballistic  vehicle  re-entering  the  earth's  atmosphere  at  high  speeds 
possesses  a  large  amount  of  kinetic  energy  which  is  converted  into  aerodynamic 
heat.  The  terrific  heat  which  surrounds  the  missile  body  requires  that  the  bgdy 
be  coated  with  some  type  of  material  that  can  withstand  this  intense  heat.  / Theo¬ 
retical  and  experimental  investigations  have  proved  that  certain  type  glasses 
makef  a  very^good  protective  covering  for  a  re-entry  vehicle. 


For  the  approximately  spherical  surface  in  the  vicinity  of  the  forward  stag¬ 
nation  point  of  an  ICBM,  this  report  presents  calculated  results  pertaining  to 
the  transient  heating  and  ablation  processes  of  26  homogeneous,  opaque,  and  non¬ 
decomposing  glass  shields,  each  shield  having  a  unique  set  of  thermal  properties. 
The  effect  of  the  thermal  properties  on  the  heating  and  ablation  results  was 
investigated  by  varying  each  thermal  property  over  a  practical  range  of  values. 

A  similar  studv  of  the  effect  of  the  thermal  properties  on  the  heating  and  abla-'*° 

*■'  -  I  n  | ,  gi  IM  I  main  I  Mil  my  - ■  . .  • 

tion  results  has  been  made  for  an  IRBM  re-entry  I  [1,  4]. 


y  fee  results  , presented  in  this  report  for  the  ICBM  re-entry  were  obtained  by 
employing  th^^mnnerica'i 7  calculation  me thodj derived  in  references  3  or  5yan3Twere 
evaluated  on  an  IBM  704  computer.  W  _ 

. '  f>  \y 


THE  TRAJECTORY  AND  BODY  CHARACTERISTICS 


The  ICBM  is  assumed  to  begin  its  re-entry  into  the  earth’s  atmosphere  at  an 
altitude  of  150  km  at  Mach  15.68  or  a  flight  speed  of  7225.9  m/sec.  Figure  2 
shows  the  flight  altitude  and  flight  speed  as  a  function  of  time  for  the  ICBM 
re-entry.  The  time  t  =  0  sec  corresponds  to  the  altitude  of  150  km.  The  re-entry 
angle  (0g)  is  114°,  measured  from  the  vertical,  and  the  angle  of  attack  (a)  that 
remains  constant  throughout  the  re-entry  is  0°.  The  horizontal  range  of  the 
missile  from  zero  time  (H  =  150  km)  to  impact  time  (H  =  0  km)  is  330.3  km.  The 
ballistic  factor  (cB)  of  the  ICBM  is  1000  lb/ft2  or  0.002  m3/kgsec2  according  to 
different  definitions. 

The  spherical  portion  of  the  glass  shield  at  the  stagnation  point  is  assumed 
to  have  a  diameter  of  0.635  m. 


MATERIAL  PROPERTIES  AND  FUNCTIONS  EMPLOYED 


Since  there  is  a  lack  of  experimental  data  at  high  temperatures  for  the 
thermal  conductivity  (k) ,  the  specific  heat  (cp),  and  the  emissivity  constant 
(e),  these  material  properties  were  treated  as  non-temperature-dependent  con¬ 
stants  for  each  calculated  glass  shield.  The  effect  of  the  variation  of  these 
material  property  constants  is,  however,  investigated  by  the  calculation  of  26 
supposedly  opaque  glass  shields  (see  Table  1)  with  each  glass  shield  having  a 
different  set  of  material  properties.  The  specific  weight  (7)  is  assumed  to  be 
2105  kg/m3  for  all  the  calculated  glass  shields  in  this  report. 

Since  experimental  data  is  available  for  the  viscosity  and  the  vapor  pressure 
as  a  function  of  the  temperature,  these  two  material  properties  of  glass  were 
treated  as  temperature-dependent  functions  for  each  calculated  glass  shield.  The 
viscosity  functions  for  Pyrex  (Corning  7900)  and  fused  silica  (quartz)  are  curve 
fits  of  curves  presented  in  Figure  9  of  reference  11.  The  vapor  pressure  function 
and  heat  of  vaporization  (hv)  for  Pyrex  were  taken  from  reference  7.  The  vapor 
pressure  function  and  heat  of  vaporization  for  fused  silica  were  taken  from 
reference  9.  When  the  constants  jj*  and  A  in  the  viscosity  and  vapor  pressure 
functions  are  equal  to  one,  the  equations  given  below  for  \i  and  pv  are  the  exact 
curve  fits  of  these  two  properties  for  Pyrex  and  fused  silica  glasses. 

The  vapor  pressure  function,  viscosity  function,  and  heat  of  vaporization 
for  Pyrex  and  fused  silica  glass  are  as  follows: 

(1)  Pyrex  glass 

(a)  the  vapor  pressure  function 

PV,P  =  10,332  Ap  exp  +  14.5  (kg/m2)  (1) 

(b)  the  viscosity  function 

/A  on  A  \  2. 0635-1 

Pp  =  0.016684  p p  exp  ( '^ofcy)  (k8  sec/m2)  (2) 

(c)  the  heat  of  vaporization 

hv?p  =  2470  (kcal/kg) 

(2)  Fused  silica  (quartz) 

(a)  the  vapor  pressure  function 

pv,F.S.  =  10>332  exp  T(°K)°  +  18'^8  (kg/m2)  (3) 
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(b)  the  viscosity  function 

oo £  Q7\2.6366- 

^F.S.  =  ^,9505  exP  t(°K)  )  (kS  sec/m2)  (4) 

(c)  the  heat  of  vaporization 
hv,F.S.  :=  3050  <kcal/k§> 

The  constants  [_ip  and  Ap  are  one  for  the  real  material  and  were  varied  in  some  of 
the  cases  studied  here. 

The  material  properties  that  were  varied  and  the  limits  within  which  each 
material  property  is  varied  in  this  report  are  as  follows: 

(1)  the  thermal  conductivity,  k  (kcal/m  °K  sec) 

0.61  x  10"4  <  k  <  48.84  x  10~4 

(2)  the  specific  heat,  cp  (kcal/kg  °K) 

0.1.0  <  cp  <  0.58 

(3)  the  emissivity  constant,  e  (-) 

0.40  <  e  <  1.0 

(4)  the  viscosity  level,  p*  (-) 

0. 10  <  Pp  <  CO 

(5)  the  vapor  pressure  level  Ap 
0  <  Ap  <  1000 

Since  the  cases  as  the  individual  material  properties  k,  cp,  and  p*  approach  zero 
and  infinity  were  discussed  in  reference  1,  they  will  not  be  repeated  in  this 
report. 

A  so-called  "standard  set"  of  material  properties  will  be  referred  to  in  this 
report,  and  they  are  as  follows: 

k  =  0.61  x  10-4  (kcal/m °K  sec) 

cp  =  0.29  (kcal/kg  °K) 
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e  =  0.8 

(-) 

p*  =  1.0 

(-) 

II 

O 

(-) 

The  "standard  set"  of  material  properties  was  employed  in  the  theoretical  study 
of  the  ablation  of  a  glass  shield  at  the  stagnation  point  of  a  re-entering  IRBM 
[1]  and  a  re-entering  satellite  [6]. 


THE  EFFECT  OF  THE  MATERIAL  PROPERTIES  ON  THE  INTERFACE  TEMPERATURE 

The  interface  temperature  T^  follows  the  same  trend  with  respect  to  time  for 
each  glass  shield  regardless  of  the  material  properties,  and  the  maximum  inter¬ 
face  temperature  occurs  at  approximately  the  same  time  for  each  shield  since 
he  »  hj_.  Therefore,  the  effect  of  the  material  properties  on  the  interface 
temperature  can  be  studied  by  the  consideration  of  the  maximum  interface  tempera¬ 
ture  (T^  max)  for  the  ICBM  re-entry.  Figure  3  shows  the  effect  of  the  material 
properties  on  the  maximum  interface  temperature.  Large  values  for  the  thermal 
conductivity  k,  specific  heat  Cp ,  emissivity  constant  e,  and  vapor  pressure  level 
A  tend  to  keep  the  surface  temperature  low;  the  viscosity  level  p*  has  very  little 
effect  on  the  surface  temperature. 

Figure  4  is  a  pictorial  view  of  a  Pyrex  glass  shield  (Case  24,  Table  1) 
showing  its  complete  temperature  and  ablation  history  for  the  ICBM  re-entry.  One 
can  see  the  movement  of  the  interface  at  the  melting  rate  -vGO (mm/sec)  toward 
the  thermally  insulated  inner  edge  of  the  glass  shield  in  this  graph. 


THE  HEAT  FLUXES 
1.  tflie  Aerodynamic  Heat  Flux 

The  compression  at  and  near  the  forward  surface  and  the  friction  between 
the  fluid  particles  as  they  flow  along  the  body  converts  the  kinetic  energy  of 
motion  into  aerodynamic  heat.  The  formulas  that  were  used  for  the  aerodynamic 
heat  flux  rate  at  the  stagnation  point  of  a  non-ablating  surface  and  the  source 
of  these  formulas  are  as  follows: 

(1)  Molecular  flow  region  [10] 


q 


aero 


=  1.054  x  10-3  p  V 

00  oo 


V2  -  1.9263  x  103  (T. 
oo  i 


T  ) 
oo 


kcal/m2sec 


(5) 
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(2)  Continuum  flow  region  for  Vqo  >  2100  m/sec  [8] 


=  74,571.6  sfcjR  (V  /u  )3'15  (~r=m  kcal/m2sec 


aero 


oo  c 


h  -  h. 


h  -  73 
e 


(6) 


(3)  For  V  <  2100  m/sec  [12] 
v  y  oo 


Vro-O^i.air  <Te  '  V  J  WTi'  <»“  ^Pr  -  1  <-715>°‘4  %  ’  V>' 

kc al /m2sec. 

(7) 

The  aerodynamic  heat  flux  for  an  ablating  surface,  e.g.,  glass  vaporization,  is 
given  by 


q  =  q  t 
aero  aero 


where  \| /  is  given  in  Reference  7  as 


*  - 


1  - 


1  -  c 

(1  -  Kx  Mni)c 


(8) 


(9) 


c 


1 


1  +  M  ( 


£e 

Pv 


(10) 


For  the  glass  shields  considered  herein,  M  =0.72,  =  0.68,  and  nj  =  0.26. 

The  relation  for  the  surface  ablation  taken  from  reference  7  is 


(7v),  = 


q  c 
^aero 


1  (he  -  h.)  [1  -  (1  -  Kx  Mni)c] 


(ID 


The  time  t  =  t*(sec)  designates  the  time  after  ablation  has  ended  when  the 
net  heat  flux 


q.  =  q  -  q  ,  +  h  (7v) .  (12) 

aero  nrad  vw  ' x 

across  the  surface  changes  its  sign  from  positive  to  negative.  Beginning  then, 
the  radiative  heat  transfer  qracj  out  of  the  shield  exceeds  the  aerodynamic  heating 
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rate  qaero»  The  time  tw  is  between  the  time  of  peak  heating  and  impact  time  and 
is  approximately  the  same  for  all  the  glass  shields.  All  of  the  individual  heat 
fluxes  for  each  glass  shield  were  integrated  over  time  t  to  t  a  tft,  since 

•JU 

t'v 

Qi(  t*)  =  f q i  dt 

o 

gives  the  maximum  amount  of  heat  that  the  shield,  in  both  the  solid  and  the  liquid 
stage,  will  absorb.  The  time  integrated  q's  to  t“  will  be  designated  Q(t'v)  in 
this  report. 


The  aerodynamic  heat  flux  Qaero^*)  that  stands  for  the  time  integrated 
aerodynamic  heat  flux  without  mass  flow  across  the  interface  is  given  for  each 
one  of  the  calculated  glass  shields  in  Table  1.  The  true  aerodynamic  heat  flux, 
Qaero(t*)»  is  given  by 


Q  (t*)  =  Q  (O  -  Q  -  + 
xaero  %aero  Bl 


h  (-/v) .  dt, 
v  v  '  1 


(13) 


where  Q^^(tw)  is  the  total  amount  of  aerodynamic  heat  blocked  up  to  time  t'v  due 
to  the  vaporization  process  of  the  glass.  Table  1  shows  that  the  total  aero¬ 
dynamic  heat,  Qaero(t*):>  changes  very  little  for  all  of  the  calculated  cases 
since  he  »  hi  during  the  period  of  significant  heating.  Thus,  the  material 
properties  of  the  glass  shield  have  very  little  effect  on  the  aerodynamic  heat 
flux  Qaero*  The  trajectory  of  the  re-entering  vehicle  principally  determines 
the  aerodynamic  heat  flux  Qaero* 


2.  The  Seat  Flux  Radiated  Away  From  the  Opaque  Glass  Shields 


The  heat  flux  that  is  radiated  away  from  the  outer  surface  of  the  sup¬ 
posedly  opaque  glass  shield  given  by 


q  ,  =  1.378  x  1CT11  €  T?  (kcal/m2  sec) 
rad  l 


(14) 


is  a  function  of  the  wall  temperature  and  emissivity  constant  e.  The  values  of 
Qrad(t*)  listed  in  Table  1  show  that  the  outer  surface  radiates  away  the  major 
portion  of  the  aerodynamic  heat  flux  Qaero(t*)  for  most  of  the  calculated  shields. 
The  radiative  efficiency  of  the  glass  shields  is  defined  as  the  fraction 
Qrad(t*)/Qaero(t*)  of  the  gross  heat  input  that  is  radiated  away  from  the  glass 
wall.  Figure  5  shows  the  radiative  efficiency  of  the  heat  protection  shields  for 
the  varied  material  properties  assumed  for  the  glass.  This  graph  shows  that  the 
radiative  efficiency  of  a  glass  heat  shield  is  optimized  by  utilizing  a  glass 
with  a  small  thermal  conductivity  and  specific  heat  and  with  a  small  vapor  pressure. 
Changes  in  the  viscosity  level  and  emissivity  constant  have  a  much  smaller  effect 
on  the  radiative  efficiency  than  the  other  pertinent  material  problems. 
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3.  The  Heat  Blocked,  (^(t"),  Due  to  the  Vaporization  Process  of  the  Glass 

The  total  amount  of  heat  blocked  due  to  the  vaporization  is  given  by 

t* 

Qb1  =  Q  -  Q  +  h  /  (-tv),  dt,  (15) 

Bl  aero  aero  v  J  '  1 

o 

where  (Qaero  ~  Qaero)  represents  the  amount  of  heat  blocked  by  the  diffusion  of 
the  vapor  across  the  boundary  layer  and 

t* 

hv  J  ("7v)i  dt 

o 

is  the  amount  of  heat  absorbed  by  the  evaporation  process  at  the  surface  of  the 
shield.  Cases  calculated  for  Ap  =  1,  3,  100,  and  1000  (Table  1),  with  all  other 
material  properties  being  the  so-called  "standard  set"  of  material  properties, 
resulted  in  46  to  88  percent  of  the  total  aerodynamic  heat  Qaero(t*)  being  blocked. 
The  amount  of  the  total  ablation  that  was  due  to  vaporization  ranged  from  53  to 
98  percent.  Reference  5  points  out  that  radiative  cooling  is  more  desirable  than 
mass  transfer  cooling  under  certain  conditions. 

4.  The  Heat  Flux  Q^(t*)  Absorbed  by  the  Glass  Shield 

The  total  amount  of  heat  that  the  molten  and  solid  layers  of  the  shield 
absorb  is  given  by 


Qi(t*>  =  \ero(tVf)  "  Wfc*>  "  " 


=  Q  (t'f)  -  Q  ,(t*)  -  h  /  (-yv) .  dt. 
^aero'  7  ^radv  7  v  \  '  l 


(16) 


This  total  heat  absorbed  Q^(t ')  per  unit  area  by  the  shield  ranges  from  3146  kcal/m2 
to  60^000  kcal/m2  for  all  cases  calculated.  The  amount  of  the  total  aerodynamic 
heat  Qaero(t*)  absorbed  by  the  shield  ranged  from  7  to  91  percent. 

Figure  6,  presenting  Q^(t*)  as  a  function  of  the  material  properties,  shows 
that  the  thermal  conductivity  and  specific  heat  affect  the  magnitude  of  absorbed 
heat  Q^(t*)  by  the  shield  much  more  than  the  other  material  properties.  The 
amount  of  heat  Q^(t'')  should  be  small  since  then  the  shield  thickness  needed  as 
a  heat  sink  is  kept  to  a  minimum.  Qj_(t^)  can  be  kept  small  by  employing  a  glass 
with  a  small  thermal  conductivity,  a  small  specific  heat,  a  high  viscosity  and 
vapor  pressure  level,  and  a  large  emissivity  at  the  surface. 
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5,  The  Amount  Qc(t*)  of  the  Absorbed  Heat  Carried  Away  by  the  Molten  Flow 

The  convective  heat  flux  Qc(t'')  is  the  amount  of  heat  that  the  flow  of 
molten  glass  carries  away  from  the  stagnation  point  vicinity.  This  heat  Qc(tx) 
carried  away  by  the  liquid  glass  is  part  of  the  heat  Qi(t*)  that  was  absorbed  by 
the  glass  shield.  Table  1  shows  that  from  0  to  79  percent  of  the  total  aero¬ 
dynamic  heat  Qaero^*)  is  carried  away  by  the  molten  flow.  The  amount  of  the 
heat  absorbed  Qi(t*)  that  is  carried  away  by  the  molten  flow  ranged  from  0  to 
96  percent.  Table  1  and  Figure  6  show  that  those  shields  which  absorbed  a  large 
percentage  of  the  total  heat  QaeroC*1*)  ^ad  most  of  this  absorbed  heat  Qj-(t*)  ; 

carried  away  by  the  molten  flow. 

The  difference  Q -^ ( )  -  Qc(t*)  is  that  amount  of  heat  the  glass  shield 
will  have  to  store  in  heat  sink  fashion  at  the  stagnation  point,  because  after 
time  t  no  more  convective  cooling  takes  place  since  melting  has  stopped.  After 
time  t*  the  heat  Qi(t*)  -  Qc(t*)  is  very  important  since  this  heat  keeps  pene¬ 
trating  into  the  shield  until  impact  time  tf.  The  fraction  R(t^v)  given  by 
(Qi(tVc)  -  Qc(t*))  /Qaero^*)  fraction  of  the  gross  heat  Qaero^*)  that  is 

stored  in  heat  sink  fashion  in  the  shield.  This  fraction  R(t^c)  varies  from  0.027 
to  0.122;  this  shows  that  the  glass  shields  have  to  store  very  little  of  the  gross 
heat . 


THE  ABLATION  OF  THE  GLASS  SHIELDS 

The  ablation  rate,  -vO0  (mm/sec),  of  the  glass  shield  is  due  to  (1)  the 
flow  of  molten  glass  and  (2)  the  evaporation  of  the  glass  at  the  surface.  The 
total  ablation  s,  i.e.,  the  layer  lost  due  to  ablation,  given  by 


for  each  glass  shield,  is  given  in  Table  1.  For  the  shields  investigated,  the 
thickness  s  of  the  heat  shield  lost  due  to  ablation  was  from  0  ,to  35  mm,  amounting 
to  less  than  1,4  inches. 

Figure  7  shows  the  total  ablation  s  as  a  function  of  the  material  properties 
of  the  glass  shield.  Figures  7,  8,  and  9  show  that  the  thermal  conductivity  k, 
and  the  specific  heat  Cp  are  the  material  properties  that  principally  determine 
the  total  ablation  s  when  the  ratio  k/7cp  is  small,  as  is  desirable.  Practical 
variations  of  the  viscosity  level  f!*,  surface  emissivity  e,  and  the  vapor  pressure 
level  A  caused  rather  small  variations  in  the  total  ablation  for  the  glass  shields 
which  had  a  small  thermal  diffusivity  k/7cp. 

The  dependency  of  s  on  the  viscosity  level  for  a  non-vaporizing  glass 
shield  is  shown  in  Figure  8  by  comparison  of  s  for  cases  8,  19,  and  20.  All  three 
cases  have  the  properties  k  =  9,7  x  10"4  kcal/m  °K  sec,  cp  =  0.29  kcal/kg  °K, 
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k/ycp  =  15.9  x  lCT7  m2/sec,  and  Ap  =  0.  The  viscosity  level  jj*p  was  varied  within 
the  limits  0.1012  <  |iXp  <  101.2  which  resulted  in  the  total  ablation  s  ranging 
between  9.3  <s  <35.7  mm.  Evidently,  s  increases  with  a  decreasing  viscosity 
level. 

The  effect  of  varying  k  and  cp  for  a  large  constant  ratio  k/ycp  is  shown 

by  comparing  Cases  21  and  22  (Fig.  8)  where  k#cp  =  40  x  10’7  m2/sec  for  both 

cases;  Case  21  has  values  of  k  and  cp  which  are  smaller  than  those  of  Case  22 

by  a  factor  of  two.  Case  21  with  the  smaller  k  and  cp  resulted  in  an  ablated 

layer  59%  larger  than  that  of  Case  22. 

Cases  2,  10,  11,  25,  and  26  in  Table  1  and  Figure  7  show  the  effect  on  the 
total  ablation  s  of  varying  the  vapor  pressure  level  Ap  between  the  limits 
0  <  Ap  <  1000,  The  total  ablation  s  increases  only  slightly  with  increasing 
vapor  pressure  level.  Table  1  shows  that  the  fraction 


of  the  total  ablation  due  to  vaporization  is  as  large  as  98%  for  Ap  =  1000,  and 
that  this  fraction  increases  together  with  the  vapor  pressure  level  A  and  vis¬ 
cosity  level  p*. 


THE  THERMAL  PENETRATION  INTO  THE  GLASS  SHIELD 
The  calculated  "semi-infinite  solution"  [3]  is  based  on  the  condition  that 
lim  T(t ,  z)  *  Tq 

Z  — >  oo 

where  T0  was  assumed  to  be  300°K  in  this  report.  For  all  the  cases  considered, 
a  thermally  insulated  inner  wall  of  the  shield  was  assumed  to  have  been  placed 
at  the  maximum  distance  (denoted  z  =  b)  from  the  final  surface  where  the  semi¬ 
infinite  solution  yielded  T  =  380°K.  For  the  givBn  re-entry,  the  time  of  maxi¬ 
mum  thermal  penetration  was  impact  time  tf.  The  true  temperature  at  the  point 
b  and  time  tf  is 

T(tf ,  b)  =  T0  +  2(380  -  T0)  =  460°K,  (18) 

which  follows  from  adding  the  temperature  profile  given  by  the  semi-infinite 
solution  to  the  reflection  of  this  profile  about  the  point  b.  This  satisfies 
the  conditions  that  &T(tf,  b)/<3z  =  0  and  that  the  temperature  of  the  thermally 
insulated  inner  wall  is  <  460°K.  The  numerical  results  show  that  this  assumption 
causes  a  negligible  temperature  increase  at  the  outer  surface. 
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The  maximum  thermal  penetration  b,  or  the  thickness  of  glass  needed  as  a 
heat  sink,  is  given  in  Table  1  for  all  of  the  calculated  cases  and  is  plotted  as 
a  function  of  the  material  properties  in  Figure  7.  The  thickness  b  varies  from 
2.6  to  21.15  mm  (Table  1).  Figure  7  shows  that  the  thermal  penetration  b  is 
predominantly  governed  by  the  thermal  conductivity  k  and  specific  heat  Cp.  The 
viscosity  level  jj"  has  some  effect  on  b,  but  not  as  much  as  k  or  Cp,  whereas  the 
emissivity  e  ,of  the  surface  and  the  vapor  pressure  level  A  had  an  almost  negli¬ 
gible  effect  on  this  thickness  b  (Fig.  7  and  8)  since  large  percentages  of  the 
total  radiation  emitted  from  the  surface  and  the  total  transport  of  material  by 
melt  flow  and  evaporation  are  concentrated  in  the  short-lasting,  high-temperature 
range  of  the  surface.  In  this  range,  small  changes  in  the  surface  temperature 
can  compensate  for  large  changes  in  the  emissivity  constant  and  in  the  assumed 
level  factors  A  and  p*  of  the  vapor  pressure  and  viscosity,  respectively,  because 
of  the  highly  nonlinear  relations'  pv  =  pv(T),  p  =  p  (T)  ,  and  qrad  =  ^radO^)  * 

Figure  8  shows  the  thermal  penetration  b  as  a  function  of  the  thermal 
diffusivity  k/K Cp  for  all  the  calculated  glass  shields;  it  is  seen  that  b 
increases  together  with  k/ycp  and  p" 

THE  NECESSARY  WEIGHT  OF  THE  GLASS  SHIELDS 

The  necessary  weight  per  unit  area  of  the  glass  shield,  defined  as  y(s  4-  b) , 
is  the  most  important  parameter  of  the  ablation  type  heat  shield  to  be  considered; 
the  sum  (s  +  b)  is  defined  as  the  necessary  thickness.  The  weight  of  the  protec¬ 
tive  glass  shield  should  be  kept  as  light  as  possible,  but  the  shield  should 
still  successfully  protect  the  re-entry  body.  The  specific  weight  y  of  the  glass 
shields  was  assumed  to  be  constant,  i.e.,  y  =  2105  kg/m3,  in  this  report; 
therefore,  the  two  constituents  that  determine  the  necessary  weight  y  (s  +  b)  are 
the  total  ablation  s  and  the  heat  sink  thickness  b. 

For  an  IRBM  re-entry  [4],  the  thermal  diffusivity  k/ycp  was  the  dominant 
material  parameter  that  determined  the  necessary  thickness  (s  +  b)  of  a  glass 
shield  in  the  vicinity  of  the  stagnation  point.  Melting  and  significant  aero¬ 
dynamic  heating  stopped  about  40  seconds  before  impact  time  for  the  IRBM  re-entry. 
Since  s  was  considerably  smaller  than  b,  the  dependency  of  (s  4-  b)  on  the  material 
properties  other  than  k/ycp  was  small. 

For  the  ICBM  re-entry  at  hand,  melting  and  significant  aerodynamic  heating 
stopped  only  15  seconds  before  impact  time;  also,  the  contribution  of  the  total 
ablation  s  to  the  necessary  thickness  (s  +  b)  is  quite  significant  except  for  the 
cases  where  the  viscosity  level  was  assumed  to  be  so  high  that  little  melting 
occurred.  Therefore,  (s  -f  b)  is  affected  by  other  parameters  in  addition  to 
k/y  Cp .  Figure  9  presents  the  necessary  weight  as  a  function  of  the  thermal  dif¬ 
fusivity  k/7Cp  for  all  the  glass  shields  calculated.  From  the  results  it  is 
concluded  that  the  following  trend  in  the  material  properties  would  provide  the 
lightest  glass  heat  protection  shield  for  the  ICBM  re-entry: 
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(1)  a  small  thermal  diffusivity  k/ycp  , 

(2)  a  high  viscosity  level,  i.e.,  a  large  resistance  to  melting, 

(3)  a  large  emissivity  €  of  the  supposedly  opaque  glass. 

The  results  show  that  the  vapor  pressure  level  A  has  a  negligible  effect  on  the 
necessary  weight  y  (s  +  b) . 


^CONCLUSIONS 


The  heat  protection  at  the  stagnation  point  of  an  ICBM  re-entering  the 
atmosphere  of  the  earth  is  investigated  in  this  report  by  using  homogeneous, 
opaque  glass  shields.  The  effect  of  the  material  properties  of  the  glass  on  the 
performance  parameters  is  analyzed  from  the  results  of  26  calculated  glass 
shields,  each  with  a  different  set  of  assumed  material  properties. 


Tie  most  important  performance  parameter  of  the  heat  protection  device  is 
the  necessary  weight  7  (s  +  b)  per  unit  area  that  is  needed  to  keep  the  thermally 
insulated  inner  edge  of  the  shield  below  an  arbitrarily  chosen  low  temperature 
of  460°K .  The  following  trends  of  the  glass  thermal  properties  would  tend  to 
minimize  the  weight  of  the  glass  heat  protection  shield  for  the  ICBM  re-entry: 


(1)  a  small  thermal  diffusivity  k/ycp, 

(2)  a  high  viscosity  level,  i.e.,  a  large  resistance  to  melting, 

(3)  a  large  emissivity  constant  e  of  the  supposedly  opaque  glass; 

this  is  not  as  important  as  the  thermal  diffusivity  or  the  viscosity 
of  the  glk'sb  . 


It  was  found  that  the  vaporization  of  the  glass  had  very  little  effect  on 
the  necessary  weight  7  (s  +  b)  .1 
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FIG.  1.  THE  REFERENCE  SYSTEM 
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FIG.  2.  THE  TRAJECTORY  OF  THE  ICBM 
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THE  THERMAL  DIFFUSIVITY  k/yc  FOR  ALL  THE  CALCULATED  GLASS  SHIELDS 
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FIG.  9.  THE  NECESSARY  WEIGHT  y(s+b)  AS  A  FUNCTION  OF  THE  THERMAL 
DIFFUSIVITY  FOR  ALL  THE  CALCULATED  GLASS  SHIELDS 
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